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A standard approach to analysis of noise-induced effects in stochastic dynamics assumes a Gaus- 
sian character of the noise term describing interaction of the analyzed system with its complex 
surroundings. An additional assumption about the existence of timescale separation between the dy- 
namics of the measured observable and the typical timescale of the noise allows external fluctuations 
to be modeled as temporally uncorrelated and therefore white. However, in many natural phenom- 
ena the assumptions concerning the abovementioned properties of "Gaussianity" and "whiteness" of 
the noise can be violated. In this context, in contrast to the spatiotemporal coupling characterizing 
general forms of non-Markovian or semi-Markovian Levy walks, so called Levy flights correspond to 
the class of Markov processes which still can be interpreted as white, but distributed according to a 
more general, infinitely divisible, stable and non-Gaussian law. Levy noise-driven non-equilibrium 
systems are known to manifest interesting physical properties and have been addressed in various 
scenarios of physical transport exhibiting a superdiffusive behavior. Here we present a brief overview 
of our recent investigations aimed to understand features of stochastic dynamics under the influence 
of Levy white noise perturbations. We find that the archetypal phenomena of noise-induced ordering 
are robust and can be detected also in systems driven by non-Gaussian, heavy-tailed fluctuations 
with infinite variance. 

PACS numbers: 05.40.Fb, 05.10.Gg, 05.40.-a, 02.50.-r, 05.40.Ca. 



I. INTRODUCTION 

Systems operating far from thermodynamic equilib- 
rium are usually subjected to the action of noise which 
may profoundly influence their performance. Over the 
past two decades the extensive theoretical and experi- 
mental studies in physics, information theory and biology 
have documented various phenomena of noise-induced 
order, noise-facilitated kinetics or noise-improved signal 
transmission and detection. Almost all research in this 
field assumes that the noise process involved can be char- 
acterized by finite variance. Yet, in many situations the 
external non-thermal noise can be described by distribu- 
tion of impulses following the heavy-tail stable law statis- 
tics of infinite variance. Notably, infinite variance does 
not exclude finite spread or dispersion of the distribution 
around its modal value. In fact, for the family of a-stable 
noises the interquantile distance can be used as a proper 
measure of the distribution width around the median. 

Dynamical description of the system and its surround- 
ings is typically carried out within a stochastic picture 
based on Langevin equations [l|, 0] • The basic equation 
of this type reads 



x(t)=f(x,t) + ((t), 



(1) 



where f(x,t) is the deterministic (and possibly time-de- 
pendent) "force" acting on the system and £(i) stands for 
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the "noise" contribution describing interaction between 
the system and its complex surrounding. If the noise can 
be considered as white and Gaussian, the above equation 
gives rise to the classical Langevin approach used in the 
analysis of Brownian motion. The whiteness of the noise 
(lack of temporal correlations) corresponds to the exis- 
tence of time-scale separation between the dynamics of 
a relevant variable of interest x(t) and the typical time 
scale of the noise. Hence white noise can be considered 
as a standard stochastic process that describes in the 
simplest fashion the effects of "fast" surroundings [!, H| . 

Although in various phenomena, the noise can be in- 
deed interpreted as white (i.e. with stationary , indepen- 
dent increments), the assumption about its Gaussianity 
can be easily violated. The examples range from the 
description of the dynamics in plasmas p|, diffusion in 
energy space, self-diffusion in micelle systems, exciton 
and charge transport in polymers under conformational 
motion and incoherent atomic radiation trapping - to the 
spectral analysis of paleoclimatic @ or economic data 0] , 
motion in optimal search strategies among randomly dis- 
tributed target sites Q, fluorophore diffusion as studied 
in photo-bleaching experiments, interstellar scintillations 
@ , ratcheting devices [13, [Hj] and many others [HJ . 

The present work overviews properties of Levy flights 
in external potentials with a focus on astonishing aspects 
of noise-induced phenomena [13, 0, EH like resonant ac- 
tivation (RA), stochastic resonance (SR) and dynamic 
hysteresis [13, E3, EM E3, H3]- In particular, the persis- 
tency of the SR occurrence is examined within a contin- 
uous and a two-state description of the generic system 
[HI composed of a test particle moving in the double 
well-potential and subject to the action of deterministic, 
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periodic perturbations and a-stable Levy type noises. In 
the same system the appearance of dynamic hysteresis 
is documented. Moreover, an archetypal escape prob- 
lem over a fluctuating barrier is analyzed revealing the 
RA phenomenon in the presence of non-equilibrium Levy 
type bath. 



II. MODEL 

The model system is described by the following (over- 
damped) Langevin equation driven by a Levy stable, 
white noise ((t) 



dx{t) 
dt 



(2) 



To examine occurrence of stochastic resonance and dy- 
namic hysteresis, the generic double-well potential with 
the periodic perturbation (see left panel of Fig. [T]) has 
been chosen 



a 2 b 4 



V(x, t) = + l x + A)£cosQi, 



(3) 



with a = 128, b = 512, An 
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FIG. 1: The generic double-well potential V(x,t) = — fx 2 + 
|a; 4 + AoxcosQt with a = 128, b = 512, A = 8 for 
t = {0,7r/2,7r} (the left panel) and the generic linear poten- 
tial slope dichotomously switching between two configurations 
characterized by differend heights H± (right panel). 



In turn, for inspection of resonant activation, the ac- 
tion of the potential V(x,t) has been approximated by 
a linear slope switching dichotomously between two dis- 
tinct configurations V±(x) — H±x (see right panel of 
Fig. [l| . The detection of RA has been carried for a rel- 
atively "high" barrier H + — 8, which under the normal 
diffusion condition guarantees the proper separation of 
time scales (the actual times of escape events and time 
of diffusive motion within the potential well (21. l22ln of 
the process. Following former studies on RA 22|. the 
lower barrier height has been chosen to take one of two 
values H- = —8 or H- = 0. A particle has been assumed 



to start its motion at the reflecting boundary x = and 
continued as long as the position fulfilled x{t) < 1. 

It should be stressed that with generally non-Gaussian 
white noise the knowledge of the boundary location alone 
cannot specify in full the corresponding boundary condi- 
tions for reflection or absorption, respectively [23]. The 
trajectories driven by non-Gaussian white noise display 
irregular, discontinuous jumps. As a consequence, the lo- 
cation of the boundary itself is not hit by the majority of 
discontinuous trajectories. This implies that regimes be- 
yond the location of the boundaries must be properly ac- 
counted for when setting up the boundary conditions. In 
particular, multiple recrossings of the boundary location 
from excursions beyond the specified state space have to 
be excluded. As it has been demonstrated elsewhere [23[ , 
incorporation of Levy flights into kinetic description re- 
quires use of nonlocal boundary conditions which, for the 
case of absorbing boundary at, say, x = 1 calls for exten- 
sion of the absorbing regime to the semifine beyond that 
point, x 1. 

The barrier fluctuations causing alternating switching 
between the high (H + ) and low (£/_) barrier configu- 
rations have been approximated by the Markovian di- 
chotomous noise [24| with the exponential autocorrela- 
tion function ([»,(*) - fa(i))]fo(f) - fa(f))]> = k( H + ~ 
H ) 2 exp(— 27 1 1 — t'\). Here 7 represents the transitions 
rates between (±) and (=p) states, and both Levy and di- 
chotomous noises have been assumed to be statistically 
independent. 

The sample trajectories x(t) have been obtained by 
a direct integration of the Langevin equation ((2J) using 
the standard techniques of integration of stochastic dif- 
ferential equations with respect to the Levy stable PDFs 

HE Q2IH HI 

«o - -Jvu^s + jasys ( 4) 

to to 
* N-l 

to 1=0 

Here NAs = t — to and Q 1S distributed according to 
the stable probability density function L Q)| g(£; a, fx = 0) 
whose corresponding Fourier transform for a ^ 1 is 



4>{k) = exp — (T Q |fc| a ^1 — «/3signfctan 
while for a = 1 this expression reads 



~2J 



4>{k) = exp 



- <j\k\ ( 1 + ij3— sign/clog | k 

7T 



(5) 



(6) 



The parameter a £ (0, 2] denotes the stability index, 
yielding the asymptotic long tail power law for the £- 
distribution, which for a < 2 is of the |^| _1 ~ Q type. The 
parameter a (a G (0, 00)) characterizes the scale whereas 
/3 ((3 S [— 1, 1]) defines an asymmetry (skewness) of the 
distribution [25l.[27l. 
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Exemplary probability density functions for symmetric 
([3 = 0) and skewed (/3 ^ 0) cases have been displayed in 
Fig.H 




FIG. 2: Symmetric stable distributions for a — 
2.0,1.5,1.0,0.5 with (3 = (from the bottom to the 
top) (left panel) and asymmetric stable densities for /? = 
1.0, 0.5, 0.0, -0.5, -1.0 with a = 1.5 (from left to right) (right 
panel) . 




FIG. 3: Dynamic hysteresis loops for various (5 with fixed a = 
1.9. The time step of the integration At = 10 -4 . Results were 
averaged over N — 10 3 realizations. Initial conditions a;(0) 
were sampled from the interval [—1.25, 1.25]. The particle 
is moving in the modulated double-well potential Q with 
a = 128, b = 512, A = 8, Q = 1, a = y/2. 



The overall behavior of the dynamic hysteresis loops 

27r/f2 

HL= J (x(t))d(cosnt) (8) 
t=o 



III. RESULTS 

A. Stochastic resonance and dynamic hysteresis 

The response of the ensemble average coordinate (x(t)) 
to external periodic perturbations (cf. Eq. ([2])) gives rise 
to the phenomenon of "dynamic hysteresis" . In order to 
register the behavior, we have performed an analysis of 
trajectories x(t) based on a two-state approximation. For 
this purpose, we have defined an occupation probability 
of being in the left (with the respect to the barrier posi- 
tion x = Xb) or right state according to the definition of 
the occupation probabilities: 

Pich(t) = Prob[a;(t) < x b ] = J p(x,t)dx = 1 -p r j g ht(*)- 

— oo 

(7) 

Note that the relative occupation of either one of those 
states can be altered not only by the modulation of the 
potential, but also by an appropriate tuning of the ad- 
ditive noise parameters. When plotted as a function of 
the periodic driving cos fit, the occupation probability 
Pright(i) exhibits a characteristic hysteresis loop. Quite 
obviously, the nonzero skewness parameter (3 introduces 
asymmetry to the dynamic hysteresis loops, cf. Fig. [3J 
As can be inferred from Fig. [31 the positive (3 "biases" 
the motion, causing the trajectories to stay more likely 
to the right of the barrier top x = 0. 



under the a-stable noises could be deduced from the in- 
spection of the trajectories of the process presented in 
Fig. For the symmetric stable noise, i.e. (3 = 0, the 
process spends, on average, the same amount of time 
in the left /right states and consequently, occupations of 
both states are equal. For non-zero and increasing |/3| a 
larger asymmetry in the distribution of residence times in 
the right /left states is registered. The asymmetry of the 
occupation probability in either one of the states is re- 
flected in the shape of the dynamic hysteresis loop, which 
with an increasing \(3\ becomes distorted into the direc- 
tion determined by the sign of the skewness parameter. 

The probability of finding the process in the right /left 
state depends strongly on the stability index a, cf. Fig. 2] 
For symmetric noises {(3 = 0) the area of the hystere- 
sis loop decreases with decreasing stability index a, see 
Fig. 2J This observation is a direct consequence of a 
heavy-tailed nature of the noise term in Eq. (p}. With 
decreasing a, larger excursions of the particle are pos- 
sible and these occasional jumps of trajectories may be 
of the order of, or even larger than the distance separat- 
ing the two minima of the potential V(x, t). A combined 
interplay of both noise parameters a, (3 may result in a 
permanent locking of the process in one of its states, see 
right panel of Fig. 2J 

In order to quantify the SR phenomenon P, [28[ we 
have used the standard measures of the signal-to-noise 
ratio SNR and the spectral power amplification r\ 
The SNR is defined as a ratio of a spectral content of the 
signal in the forced system to the spectral content of the 
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noise ((t): 

SNR = 2 



lim 



S(u>)du> 



/S N (Q). (9) 



Here Jq^^ U S(u>)duj represents the power carried by the 
signal, while Sjv(fi) estimates the background noise level. 
In turn, the spectral power amplification (rf) is given by 
the ratio of the power of the driven oscillations to that 
of the driving signal at the driving frequency f2. Closer 
examination of these quantifiers reveals that rj and SNR 
behave in a typical way both in the continuous model, 
cf. Fig. [6j and in its two-state analogue (results not 
shown). The spectral amplification r\ has a character- 
istic bell-shape form indicating detection of stochastic 
resonance within the interval of suitably chosen noise in- 
tensity a 2 . For example, if the stability index is set to 
a = 1.9 and trajectories are simulated with extremely 
weak noise intensities a 2 , the periodic signal is not well 
separated from the noisy background and consequently 
SNR stays negative. 

However, SNR becomes positive with increasing noise 
intensity what indicates emerging separation of the signal 
from the noisy background. In Fig. [S] results for rj and 
SNR analysis (in a continuous SR model) with various 
stability index a are presented. Note, that for a given 
a the power spectra for ±/3 are the same. Therefore, 
stochastic resonance quantifiers derived from the power 
spectra are equivalent for ±/3 with the same value of a. 
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FIG. 5: Trajectories of the a-stable process for the two-state 
and the continuous model constructed by use of Eq. <(2j) . The 
time step of the integration At = 10 -5 , the scale parameter 
a — \/2. The frequency of the cycling voltage Q, — l. Noise 
parameters as indicated in figures. 



Fig. [5J At decreasing values of the stability index a, 
the maximum of rj(a 2 ) drops and shifts towards higher 
values of the noise intensity. Obviously the decrease of 
a weakens the stochastic resonance and reduces system 
performance. The diminishmcnt of spectral amplification 
for a < 2 indicates that the input signal is worse repro- 
duced in the recorded output. This behavior is easier 
detectable for symmetric noises (cf. left panel of Fig. [5]) 
than for asymmetric ones (see the right panel) and can 
be readily explained by the analysis of exemplary trajec- 
tories, see Fig. [5l Sharp spikes clearly visible in the right 
corner panel of Fig. O are due to the heavy-tailed nature 
of the Levy stable distribution and they are more pro- 
nounced the smaller the index a becomes. Their presence 
indicates that for a sufficiently small a the trajectory be- 
comes discontinuous and, on average, switches between 
left and right wells of the potential V(x,t) are realized 
by sudden long-jump escape events fairly independent of 
the periodic driving. Although the SNR is less sensitive 
to variations in noise parameters (cf. Fig. 0, the shape 
of this function flattens for decreasing values of a thus 
hampering detection of the resonant value of the noise 
optimal intensity for which the SR phenomenon is most 
likely perceived. 



FIG. 4: Dynamic hysteresis loops for various a with /3 = 
(the left panel) and f3 — — 1 (the right panel). Parameters of 
the simulation like in Fig. [3] 

Switching from the continuous model to the two-state 
approximation (cf. Fig. 5) does not change the qualitative 
behavior of the SR quantifiers. [32| However, the slope 
of the decaying part of the signal-to-noise ratio is more 
flat for the two-state description than for the continuous 
model (results not shown). 

Our numerical analysis implies that rj is more sensi- 
tive to the variation of noise parameters than SNR, see 



B. Resonant Activation 

As already discussed in the preceding sections, the 
stochastic kinetics driven by additive non-Gaussian sta- 
ble white noises is very different from the Gaussian case. 
For a < 2, a test particle moving in the linear poten- 
tial can change its position via extremely long, jump 
like excursions. This in turn requires the use of nonlocal 
boundary conditions in evaluation of the mean first pas- 
sage time (MFPT) [H [H, H3, HH . In this paragraph the 
above issue is taken care of when generating first passage 
times (FPTs) by Monte Carlo simulations. Simulated 
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FIG. 6: SPA in arbitrary units (the top panel) and SNR (the 
bottom panel) for the perturbed generic double- well potential 
driven by additive Levy noises with (3 = (the left panel), 
(3 = —1.0 (the right panel) and various a. The simulation 
parameters like in Fig. [5] Lines are drawn to guide the eye. 



trajectories are representative for a motion of a test par- 
ticle over the interval [0, 1] (sec right panel of Fig. [T]) influ- 
enced by independent dichotomous switching of the po- 
tential slope and subjected to additive white Levy noise. 
A particle starts its motion at x = where the reflecting 
boundary is located. The absorbing boundary is located 
at x — 1 meaning, that the whole semi-axis [1, oo) is as- 
sumed absorbing, yielding zero PDF p(x, t) = for all 
x > 1. 

From the ensembles of collected first passage times we 
have evaluated the mean values of the distributions and 
our findings for MFPT are presented in Fig. [7] Integra- 
tion of Eq. |(2} was performed for a series of descending 
time steps of integration At = {1(T 2 , 1(T 3 , lfT 4 , 1CT 5 } 
to ensure that the results are self consistent. 

The left panel of Fig. [7] displays behavior of derived 
MFPTs as a function of the stability index a and the rate 
of the barrier modulation 7. In the right panel sample 
cross-sections of the surfaces MFPT (a, 7) are drawn. For 
a = 1, or adequately a ~ 1, the procedure of simulating 
skewed ((3 ^ 0) stable random variables becomes unsta- 
ble. This can be well explained by examining the form of 
the additive noise-term characteristic function cf>(t) - see 
Eqs. ([5]) and ([6]). The exponential functions are no longer 
continuous functions of the parameters and exhibit dis- 
continuities when a — 1, (3 7^ 0. Therefore, for the clarity 
of presentation, these parameter sets have been excluded 
from consideration. 

The depicted results indicate appearance of the RA 
phenomenon which is best visible when the potential 
barrier is switching between two barrier configurations 
characterized by H± ~ ±8. However, the dependence of 
RA on noise parameters is highly nontrivial. The over- 
all tendency in kinetics is similar to cases described in 
former paragraphs: heavier tails (a < 2) in distribution 
of noise increments C(t) result in stronger discontinuities 




I°8ioT 

FIG. 7: MFPT(a, 7 ) for H ± = ±8 with (3 = l,a = 1 (the 
left panel) and sample cross-sections MFPT(7) for various a: 
'+': a = 0.2; 'x': a = 0.8; '*': a = 0.9; a = 1.1 and 
a — 1.7 (the right panel). The results were calculated by 
direct integration of Eq. ([2]) with the time step At = 10 -5 and 
averaged over N = 10 3 realizations. Horizontal lines represent 
asymptotic values of MFPT(7 -> 0) and MFPT(7 -> 00). 
They have been evaluated by use of the Monte Carlo method 
with At = 10~ 5 and averaged over N = 5 x 10 3 realizations. 
Error bars represent standard deviation of the mean. 



of trajectories causing the RA effect to become inaudible 
and fading gradually. Notably, for totally skewed addi- 
tive noise (3=1 which acts in favor of the motion to 
the right, the RA seems to disappear for a = 0.9. It 
reappears again for smaller values of a, when the driving 
Levy white noise becomes a one sided Levy process with 
strictly positive increments which tend to push trajecto- 
ries to the right. 

To further examine the character and distribution of 
escape events from x = 0, we have analyzed survival 
probabilities G(7, t) constructed from generated trajec- 
tories at fixed values of frequencies 7. For better com- 
parison with the Gaussian RA scenario, the summary of 
results is displayed in Fig. [5] with sets of data relating to 
a = 2 and a = 0.9. The left and right panels of Fig. [8] 
present behavior of MFPT as the function of the barrier 
modulation parameter 7, the survival probability G(j,t) 
and exemplary cross-sections of the survival probability 
surface. The insets depict behavior of the survival prob- 
ability G(7,t) at short time scales (small t). 

At low rates of the barrier switching process (small 7), 
two distinct time scales of the barrier crossing events can 
be observed. The fast time scale corresponds to passages 
over the barrier in its lower state, while the large time 
scale is pertinent to the slower process, i.e. the passages 
over the barrier in its higher energetic state. This effect 
is well pronounced after the RA phenomenon sets up. 
The presence of the two time scales for small 7 and only 
one time scale for large 7 explains the asymptotic behav- 
ior of MFPTs. Namely, for low values of the switching 
rate 7, MFPT is an average value of MFPTS over both 
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FIG. 8: Comparison of the behavior of MFPTs for a = 2 
(left panel) and a = 0.9 (right panel) for H± = ±8. MFPTs' 
curves (upper panel), survival probability surfaces G(7,i) 
(middle panel) and sample cross-sections of G(7, i) surface 
(lower panel) for small ('+'), resonant ('x') and large ('*') 
7-values are depicted. Simulations parameters as in Fig. [7] 
Note the log scale on z-axis (middle panel) and y-axis (lower 
panel) . 



configuration of the barrier. With increasing frequency 
7 the distinguishable time scales coalesce and disappear. 
This is due to the fact that the barrier changes its height 
multiple times during the particle's motion so that result- 
ing MFPT describes kinetics over the average potential 
barrier. 



IV. DISCUSSION 

We have examined the influence of various types of 
stochastic a-stable drivings on dynamic properties of a 
generic two-state model system. Although our primary 
interest was to understand the effects of Levy-type driv- 
ings with the stability index < a < 2, the numerical 
analysis as applied in these studies remains valid for any 
set of parameters characterizing a-stable excitations. 

The response to external or/and parametric perturba- 
tions has been examined by analyzing qualitative changes 



in system's dynamic behavior as expressed in the onset 
of resonant activation, stochastic resonance and dynamic 
hysteresis. Due to the inherent symmetry of the poten- 
tial, the SR quantifiers constructed at a given value of 
the stability index a have been the same for ±/3. On the 
other hand, the asymmetry of the driving noise has been 
shown to influence strongly the population of states and 
therefore has affected mostly the appearance and perfor- 
mance of the dynamic hysteresis. 

The system efficiency has been described by standard 
SR measures, i.e. signal-to-noise ratio (SNR) and spec- 
tral power amplification (77). Those quantifiers have been 
shown to behave in a typical way. In other words, by tun- 
ing the noise intensity cr 2 , the maximum in the SNR(er 2 ) 
and rj(a 2 ) could be detected, thus proving that stochas- 
tic resonance is a robust phenomenon which may be ob- 
served also in systems subjected to the action of impul- 
sive jump-noise stochastic processes. Decrease in stabil- 
ity index a results in larger jump-like excursions of the 
process x(t) and consequently causes weakening of SR. 
Out of two SR measures, the spectral power amplifica- 
tion has turned out to be more sensitive to variations of a 
than SNR. A more pronounced drop in system efficiency 
has been observed for symmetric noises with (3 = 0. 

The nonequilibrated, non-thermal Levy white noise af- 
fects also a paradigm scenario of escape kinetics. By nu- 
merically implementing the set up boundary conditions 
for the problem, we have investigated the statistics of 
escape times over the dichotomously fluctuating barrier. 
The manifestation of the RA phenomenon has been an- 
alyzed within a certain frequency 7 regime pointing that 
by a continuous readjustment of the external noise para- 
maters the resonant activation can be either suppressed 
or re-induced. 

Levy white noise with a ^ 2 extends a standard Brow- 
nian noise to a vast family of impulsive jump- like stochas- 
tic processes. Our studies document that dynamical sys- 
tems driven by such sources can also benefit and display 
a noise-enhanced order. 
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